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We study exclusive B 3 decays to final states with n and r{ , induced by the rare b -> sl+r and 
b — * svv transitions. Differential decay rates and total branching fractions are predicted in the 
Standard Model, adopting the flavour scheme for the description of the r\-r] mixing. We discuss 
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I. INTRODUCTION 

The study of B s properties and decays plays an impor- 
tant role in the exploration of the Standard Model (SM) 
and in the searches for new Physics phenomena, and a 
great experimental effort is devoted at present and fore- 
seen in the near future at the hadron colliders and at the 
B factories running at the T(5S) peak [1 . In particular, 
B® — B® oscillations provide complementary information 
with respect to K° and Bd systems for the analysis of 
CP violation. CDF and DO Collaboration at the Fer- 
milab TeVatron have recently carried out measurements 
of the B® — B® mixing phase, by an angular analysis of 
the final state J/ip (f>. The measured phase is larger than 
the SM prediction (although with a sizeable error) [2J. 
If confirmed, this result would represent an evidence of 
Physics beyond SM. For this reason, new measurements 
are foreseen, involving also the final states J/t/> /o(980), 
J/ijj 77 and J/ip 77'; the branching ratio of B s — > J/iprj has 
been recently measured by the Belle Collaboration using 
the 77—5-77 and 77 — > 7r + 7r°7r~ modes to reconstruct 77 
mesons [3]. 

B s is of prime interest also for several rare decay 
modes, namely those induced by the b — > s transition, 
that are potentially important for detecting new Physics 
effects. Here we focus on the decays into 77 and 77' and 
a pair of leptons, either £ + i" (with £ = e, n and t) or 
vv. Our aim is to give predictions for several observables 
in the Standard Model, discussing the theoretical uncer- 
tainties related to the 77 — 77' mixing and to the hadronic 
matrix elements in the decay amplitudes. 

A second purpose is to consider a specific new Physics 
scenario and study how the various observables devi- 
ate from SM. The chosen framework is the Appelquist- 
Cheng-Dobrescu (ACD) model with a single universal ex- 
tra dimension (UED) j?]. The model is a minimal exten- 
sion of SM in 4 + 1 dimensions, with the extra dimension 
compactified to the orbifold S 1 /Z2 and the fifth coor- 
dinate y running from to 2ttR, y — and y = ttR 



being fixed points of the orbifold. The fields are allowed 
to propagate in all 4 + 1 dimensions, hence the model 
belongs to the class of universal extra dimension scenar- 
ios. One of its motivations is the possibility of naturally 
providing candidates for the dark matter, an issue of fun- 
damental importance . 

In the ACD model the SM particles correspond to 
the zero modes of fields propagating in the compacti- 
fied extra dimension. In addition to the zero modes, 
towers of Kaluza-Klein (KK) excitations are predicted 
to exist, corresponding to the higher modes of the fields 
in the extra dimension; such fields are imposed to be 
even under a parity transformation in the fifth coordi- 
nate P5 : y — > — y. On the other hand, fields which are 
odd under P5 propagate in the extra dimension with- 
out zero modes, and correspond to particles without SM 
partners. 

The masses of KK particles depend on the radius R 
of the compactified extra dimension, the new parameter 
with respect to SM 42J. For example, the masses of the 
KK bosonic modes are given by: 

n 2 

m l = m o+-jp 77,= 1,2,... (1.1) 

mo being the mass of the zero mode, so that for small 
values of R, i.e. at large compactification scales, the KK 
particles decouple from the low energy regime. Another 
property of the ACD model is the conservation of the 
KK parity (— l)- 7 , j being the KK number. KK parity 
conservation implies the absence of tree level contribu- 
tions of Kaluza Klein states to processes taking place at 
low energy, /j <C 1/R, forbidding the production of a sin- 
gle KK particle off the interaction of standard particles. 
This permits to use precise electroweak measurements 
to provide a lower bound to the compactification scale: 
1/R > 250 - 300 GeV [5]. Moreover, this suggests the 
possibility that the lightest KK particles are among the 
dark matter components, namely the n — I Kaluza-Klein 
excitations of the photon and neutrinos [Sj [7] . 
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Since KK modes can affect the loop-induced processes, 
Flavour Changing Neutral Current (FCNC) transitions 
are particularly suitable for constraining this new Physics 
scenario, and indeed many observables are sensitive to 
the compactification radius in case, e.g., of processes 
with B and A fc EJ [10j [HJ [12] . Here we consider rare 
B s decays into rj and 77' mesons: B s — > rj^£ + £^ and 
B s — » rj^vv, described by the effective Hamiltonian re- 



ported in Section |Hj We discuss in Section III the role 
of the 77-77' mixing and of the B s — > 77^' form factors, 
and present predictions for the various modes in the fol- 
lowing two Sections, with comments on the feasibility of 
measuring properties of these processes. Before conclud- 
ing, we discuss in the ACD model the 1/R dependence of 
the bs Unitarity Triangle, i.e. the condition among the 
Cabibbo-Kobayashi-Maskawa (CKM) elements involved 
in B s decays, and in particular the possible value of the 
CP violating phase (3 S in this model. 



II. EFFECTIVE HAMILTONIAN FOR 

B a -f ri^e+r AND B s 



(1) - 
r] y 'vv 



In the Standard Model, the effective AB = -1, AS = 
1 Hamiltonian describing the transition b — > s£ + £~ can 
be expressed in terms of a set of local operators: 



b^s£+i- 



= -4 



r 10 



(2.1) 



where Gf is the Fermi constant and are elements of 
the CKM mixing matrix (terms proportional to V u f,V* s 

v ub v: s 1 



are neglected since the ratio 



v tb v t * 



is O(10" 2 )). The 



operators Oi are written in terms of quark and gluon 
fields: 
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(2.2) 



with a, (3 colour indices, 



1±75 



6, and 



and Og denote the electromagnetic and the gluonic field 
strength tensor. 0\ and O2 are current-current oper- 
ators, 3 ,...,0 6 QCD penguin operators, 7 and O s 
magnetic penguin operators, O9 and O10 semilcptonic 
electroweak penguin operators. The Wilson coefficients 
in (2.1l have been computed at NNLO in the Standard 
Model [13J. The operators 0\ and O2 contribute to the 
the final state with a lepton pair through a cc contribu- 
tion that can give rise to charmonium resonances J/ip, 
i/j(2S), etc. The resonant term can be controlled and 
subtracted by appropriate kinematical cuts around the 
resonance masses. Since the Wilson coefficients C3 — Cq 
are small, the contribution of only the operators O7, Og 
and O\o can be kept for the description of the b — > sl + £~ 
transition, with a modification of the Wilson coefficient 
C7 described below. 

The SM effective Hamiltonian for b — > svv: 

G 

-V tb V t * s r lx X{x t )0 L = C L L 



Hb^s 



V2 2Trsm 2 (6 w ) 



involves the operator 



Or 



s y (1 _ 75 )fa/ 7jll (i _ 75 )jy 



(2.3) 



(2.4) 



Qw is the Weinberg angle; the function X{xt) (xt 



with m t the top quark mass) has been computed in |14j 
and [T5J [TB], while the QCD factor 77^ is close to one 
[T51 [TBI [T7] , so that we use 77^ = 1. 

In t he A CD m odel no operators other than those in 
(pHi ([2~2l and d2~3l contribute to b -» s£+£~ and 



b — > svv. The model belongs to the class of Minimal 
Flavour Violating (MFV) models, and the effects beyond 
SM are only encoded in the Wilson coefficients of the 
effective Hamiltonian [5]. KK excitations modify the co- 
efficients Ci and Cl, which acquire a depence on the 
compactification scale 1/R. For large values of 1/R, due 
to decoupling of massive KK states, the coefficients Ci 
and Cl reproduce the Standard Model values and the 
SM phenomenology is recovered. 

The Wilson coefficients can be expressed as functions 
F(xt, 1/R) generalizing the SM analogues Fo(xt): 



F(x t ,l/R) = F {x t ) + J2F n (x t ,x n ) , 



(2.5) 



n=l 



With Xn — 



R 2 M"L_ 



Remarkably, the sum over the KK 



— [7**, 7"]; e and g s are the electromagnetic and the strong 
coupling constant, respectively, and F^ v and G a in O7 (|2.3[) 



contributions in ( |2.5[ ) is finite, a consequence of a gener- 
alized GIM mechanism [5]; the SM results are recovered 
for R — > 0, since F(a;t, 1/R) —> Fo(xt) in that limit. 

For 1/R of the order of a few hundreds of GeV the 
coefficients differ from the Standard Model values, and 
the physical observables are predicted to be different than 
in SM. For the exclusive decays, however, it is important 
to study if this effect can be picked up, or it is obscured 
by other uncertainties, in particular those affecting the 
B s and 77, rj' matrix elements of the operators in (2.11, 
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III. 77 MIXING AND B 3 
FACTORS 



FORM 



In the determination of the B s — » rf 1 ' matrix elements 
we need to account for the rj — rf mixing. This is usu- 
ally described in two schemes: the singlet-octet (SO) and 
the quark flavour (QF) basis; in each scheme two mixing 
angles are involved [15] . 

In the SO basis one defines the rj and ^'-vacuum matrix 
elements of axial-vector currents: 



J a jP{p))=iftp„ (a =1,1 



J f/r, 

where P = 77, rf and J^ 5 , J® 5 are SU(3)f singlet and 
octet axial- vector quark currents. The four hadronic pa- 
rameters fp in (3.1 1 can be written in terms of two an- 
gles: 9 1 and 9g, and of two decay constants f\ and fg 
of a pure singlet and octet flavour state. In this scheme, 
non- vanishing values of the angles #1,8, as well as the 
difference f$ ^ /„. are SU(3)p breaking effects, and the 
contribution of the axial C/(l) anomaly is encoded in the 
decay constant f\. 

On the other hand, in the QF basis one defines the 
axial-vector currents: 



(3.1) 



and the matrix elements: 



^5 = ^7 5 s (3.2) 



(3.3) 



Also in this case the four parameters fp can be expressed 
in terms of two angles ip q and <p s , and of two decay con- 
stants fg and f s of states without and with strangeness, 
respectively |18j . The difference between the mixing an- 
gles (f q — (p s is due to OZI-violating effects and is found 
to be small (<p q — ip s < 5°), so that it has been proposed 
that the approximation of describing the r] — rf mixing in 
the QF basis and a single mixing angle is convenient |18j . 
The simplification p q ~ ip s ~ ip is supported by a QCD 
sum rule analysis of the decays 4> — > 777 and <fi — ► v'l [S3 ■ 
Here we adopt the quark flavour basis and define: 



1 

71 



I dd)) 



(3.4) 



so that the 77-77' system can be described in terms of the 
mixing angle ip = ip q = (p s : 



COS ip \TJ q ) 

sin ip\r} q ) 



- SlSHp\T] s ) 

■cos<p\r} a ) 



(3.5) 



There are several ways to measure tp, namely through 
the radiative transitions involving a light vector me- 
son V, such as V — > 77(^7 or 77' — > V7, or studying 
the two-photon decays rf'' — > 77 [20 . A precise re- 
sult has recently been obtained by the KLOE Collabora- 

. . T(6 -> 77' 7 ) 

tion which, measuring the ratio 



and assum- 



F(0 777) 

ing the flavour basis with a single mixing angle, quotes: 



ip = (41.5 ± 0.3 stQt ± 0.7 syst ± 0.6 t fc)° EQ- We use this 
value in our study. In the same analysis, KLOE has also 
allowed for a g lue content in 77', modifying the second 
relation in (3.51: 



1 77 ) = cosipGSmip\r) q ) 
+ sin ipQ \gluons) 



cos ipc cos ip \T] S 



(3.6) 



where tp G is the mixing angle for the glue contribution. 
In this case, considering also the results for the decay 
widths r(?7' — * 77), r(?7' — > fr/) and r(r/ — > W7), to- 
gether with r(vr° 77) and T(u tt°7), KLOE ob- 
tains: p = (39.7 ± 0.7)° and cos 2 <^ G = 0.86 ± 0.04. We 
comment below on how this measurement affects our pre- 
dictions. 

The mixing parameters are useful to determine the 
other quantities needed for the description of B s — > 
rf'} l + 1~ {uv) transitions, the B s and rf'^ matrix elements 
of the operators in (2.1 1, (2.3 1. As usual, such matrix el- 



ements are parameterized in terms of form factors: 
< rf^{p')\s lfl b\B s (p) >= (p + p%*f '^'V) 



M 2 B -M 2 



(3.7) 



( q = p-p>,F 1 B ^ i '\0)=F< 



fl '^ W (Q)) and 



< rf'\p')\s 1 a^q»b\B s (p) >= 

(p+j/)^-(M|.--a^w)?M 



F. 



(3.8) 



No QCD sum rule or lattice QCD calculations of B s — > 
77^' form factors (which probe the ss content of 77 and 77') 
are available, yet. The quark flavour scheme allows to 
relate the B s — > rj s form factors to the B — > K ones, so 
that: F B ^ = -simpF B ^ K and F 5 ^' = cosipF B ^ K 
(for a generic form factor F), keeping the physical masses 
of B s , 77 and 77'. It is possible to estimate the uncertainty 
connected to flavour symmetry breaking, considering re- 
lations holding in the heavy quark limit and in the chiral 
limit for heavy-to-light transition form factors [22] . Tak- 
ing as an example the form factor F\, one can write for 



In 



2f P (E P + A H ) 



(3.9) 



where H and P are a heavy and a light pseudoscalar 
meson, respectively, and Ep is the energy of the P me- 
son in the H rest frame. F is the heavy meson decay 
constant in the heavy quark limit, which is, at leading 
order in the heavy quark mass toq, independent of it 
(modulo logarithms), g describes the effective coupling 
H*HP, H* being the vector meson with the same quark 
content as H; Ah — Mh* — Mjj- If the main flavour 
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breaking effects are in Mh, Ah, Ep and fp, the ratio 
pB.-^n, /p^ K , using f K = 159.8 ± 1.4 ±0.44 MeV g3] 
and / s = (1.34 ± 0.06)/^ [18], differs from 1 by about 
10%, a reasonable estimate of SU(3)p breaking in the 
form factors. 

In the following, we use two different sets of form fac- 
tors, both obtained by QCD sum rules [24]. We re- 
fer to the first set, computed using short-distance QCD 
sum rules, as set A [25], and to the second set, com- 
puted by light-cone QCD sum rules, as set B [26 . The 
comparison allows a discussion of the uncertainty re- 
lated to the hadronic matrix elements. In both sets 
the error of the form factors is given at q 2 — and 
then extended to the full range of momentum transfer: 
Am 2 < q 2 < (M Bs -M v (,)) 2 . 

We also consider a third determination obtained apply- 
ing light-cone QCD sum rules within the Soft Collinear 
Effective Theory (SCET) [57]. It is interesting to con- 
sider this framework, since it allows to express the form 
factors in terms of a single universal function £, Ba ^ p (a 
light flavour index) through the relations, holding in the 
heavy-quark limit and in the large energy limit of the 
light meson P: 



F 



F 1 



F- 



* P (<z 2 ) 
* P (<z 2 ) 

<P (<z 2 ) 



2E' 



M 



e^ F {n+-p') (3.10) 
(n + -p') . 



In (3.101 n + is a light-like four- vector with components 
ri + = (1, 0, 0, —1), and p' is the momentum of the P me- 
son, with E' = p'° . At first order in the light pseudoscalar 
meson mass the following relations hold: 



E' 



2E' 



Mp 
2M Ba 



(3.11) 



We refer to this third parameterization as to set C. Devel- 
oping light-cone QCD sum rules in the framework of the 
SCET, in Ref.[28] the universal function ^ B ^^ relative to 
B — > 7r form factors was derived. By the same method, 
a determination can be obtained in the case of kaon, and 
the resulting function ^ B ^ K can be parameterized as: 



e- K {n + -p')=e- K { mB ) 

with 



^(m B )= 0.335™ 



a = 2.418 



094 

b = 13.765 



7 a + c ( n + -P') 

(3.12) 



c= 0.154. (3.13) 



The B s — > -q^ form factors belonging to set C can be 
written in terms of the universal function: 



Ft 



sin tp £ 



F 



IE' 



- sin tp 



■B^K 



F 



B a ^ V 



M 
M 



- sin tp 



M, 



n 



cos tp £ 



B^K 



(3.14) 



Ft 



B,->ri' 



IE' 



cos tp 



■B^K 



F 



B,->ri' 



M f 



M, 



cos tp ■ 



V' C B^K 



M f 



keeping the physical masses of the particles involved in 
the transition. Although such relations are valid in the 
heavy quark limit at small values of the momentum trans- 
fer q 2 , we extrapolate them in the full kinematical range 
of B s — > rf 1 ' transitions. 



Relating B s 



V 



W form factors to the B — > K ones 



in the QF scheme implies the neglect of flavour-singlet 
contributions, for example the annihilation contribution 
involving the strange quark spectator in B s and produc- 
ing rf'^ through two gluons. This contribution has been 
considered as possibly responsible of the anomalous pat- 
tern of the branching fractions of nonleptonic B and D 
to ?/ decays. In particular, one could expect that such 
a kind of contributions are more important for rj' than 
for 77, due to the coupling of 77' to gluons driven by the 
U(1)a anomaly in QCD (at the origin of the large 77' mass 
compared to the masses of light pseudoscalar mesons) . 

However, the size of this contribution is not firmly es- 
tablished. Denoting this gluonic contribution as Fq and 
the quark contribution as F q , a generic form factor F can 
be parameterized as: 



F 



B,-H) _ 



pBs^rf 



7-1 fl 

--hq—smtp- 
Jk 

F q ^-costp- 

JK 



Fr, 



Fa 



1_U_ 
3 fK 
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cos tp 



sin tp - 



Vsf K 



sin tp) 



cos tp) 



(3.15) 



in terms of the decay constants f q = (1.02±0.02)/ T [TB] . 
f s and fx', eqs.(3.15l show that the main effect is for the 



rj form factors. In case of B meson transitions, various 
analyses conclude that the singlet term Fq could be size- 
able, thus possibly contributing to the processes involving 
7/ [25] EH EH], however the uncertainties are large. To 
estimate such contributions the matrix elements: 



<0\A a (z)Ai(0)\v { 'Hp)> 



(3.16) 



are needed (AjJ is the gluon field, and a gauge factor has 
been dropped). These matrix elements can be written 
in terms of distribution amplitudes (DA) expanded in 
Gegenbauer polynomials. However, already the first co- 
efficient (denoted by B 2 ) of the expansion of the twist- two 
distribution amplitude (fP (assumed to be the same for 77 
and T]') is uncertain: for example, £?2(l-4GeV) = 4.6±2.5 
results from a combined analysis of the 7/777 form factor 
and of the inclusive decay Y(IS) — > rjX [35]. In cor- 
respondence to this value a contribution of about 5% is 
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estimated for the form factor F 1 B ^ r ' (0) [31], and only if 
B 2 is varied in a wider range (B 2 = 0±20) up to 0(20%) 
corrections are found in case of rf , while the corrections 
remain small in the case of r\ meson. In the framework of 
SCET these contributions arc found to be consistent with 
zero |33j . A constraint of the singlet contribution using 
the widths of the semileptonic decays B + — > rj(. + vg and 
B + — > rf IrVf, is hampered by the present experimental 
errors. 

One can investigate the presence of gluonic contri- 
butions to the form factors governing the semileptonic 
modes D s — > r]t^vi and D s — » rft^vi. The mea- 
surements BR(D S -> 7){+v t ) = (3.2 ± 0.5) x 10~ 2 and 
BR(D S -> rf£ + v e ) = (1.12±0.35) x 10~ 2 [23 can be com- 

BR(D a ^>n'Z + ve) 



bined in the ratio: Rd, 



0.35 ±0.12. 



BR{D,^ql+v t ) 

This value for R Bs is reproduced by central values of the 
mixing angle: ip = 41.5° and of the ratio of the two pieces 



contributing to Ff 3 ~ >11 written in (3.151. jrjffj 
thus pointing to small values of the singlet terms. 

A possibility for a phenomcnological analysis of B a — > 
T], rf transitions could consists in considering the effects 
of an arbitrarily added flavour singlet contribution to the 
form factors |34j . We prefer to provide predictions which 
do not include additional terms in the form factors. The 
predictions are accurate in case of ij; deviations in the 
modes with rf would first prompt investigations focused 
on the singlet contributions. 



f g (o) _ 



0. 



IV. B s -> nt+r AND B a 



The observables in these modes are the differential 
and total decay rates. The expression for the differen- 
tial width: 



dT(B s -> rfH+tr) _ G 2 F \V tb V*\ 2 a 2 A 1/2 (M 2 s , M 2 ()) ,g 2 ) / 4m 2 1 
dq 2 ~ 2 9 tt 5 M| s V q 2 3<j 



2 p(q 2 ) 



(A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz) involves the function 

p(q 2 ) = 6m 2 (M 2 Bs - M 2 (l) ) 2 \b(q 2 )\ 2 + X(M 2 Bs , M 2 (0 , q 2 ) [(2m 2 + q 2 )\c(q 2 )\ 2 - (4m 2 - q 2 )\~a{q 2 ){ 



(4.1) 



(4.2) 



where a, b and c are combinations of form factors and 
Wilson coefficients: 



~a(q 2 ) 
b(q 2 ) 



C 10 F^'(q 2 ) 

„C) 



Cio F 



C,F^-{q 2 ) 
2(m b + m s )Cj ff£ ^- 



(4.3) 



F^ ( "(q 2 ) 



M Bs + M v u 



In the expression of c in (4.3) we introduced the coeffi- 
cient , which is a renormalization scheme indepen- 
dent combination of C^,Cs and C 2 , given by a formula 
that can be found, e.g., in [§]. Using: rrib — 4.8 ± 0.2 
GeV, m s = IO4I34 MeV, together with the values: 
M B o = 5366.3 ± 0.6 MeV, t b o = (1.470±g;g|?) x 10~ 12 
s, V tb = 0.999 and V ts = (38.7 ± 2.3) x 10~ 3 quoted by 
the Particle Data Group [53], we obtain the following 
branching ratios in the SM, depending on the set of form 
factors: 



BR{B S -> r]£ + e-) 



(1.2 ±0.3) x 10 
(2.6 ±0.7) x 10 
(3.4 ± 1.8) x 10 



set A 
setB 
setC 



(4.4) 



BR{B S -» rftr) = 
for I = e, fx, and 



(1.1 ±0.3) x 10~ 7 set A 
(2.2 ±0.6) x 10~ 7 setB 
(2.8 ± 1.5) x 10~ 7 setC 



BR{B S 



BR(B S 



f]T + T 



jfr + T 



(3 ±0.5) x 10"* 
(8 ±1.5) x 10-* 
(10 ±5.5) x 10- 



(4.5) 



set A 
setB 
setC 



(4.6) 



(1.55 ± 0.3) x 10" 8 set A 
(3.85 ±0.75) x 10" 8 setB 
(4.7 ±2.5) x 10~ 8 setC 



The results obtained using the form factor £ 



b^k - 



(4.7) 
in the 

SCET framework are affected by the largest uncertainty 
and encompass the predictions based on the set A and 
set B. We consider these results as the most conservative 
determinations of the decay rates, as well as for the other 
observables. 

The results from set A and B are compatible within the 
errors, with lower central values for set A. The uncertain- 
ties, estimated at the level of 30% in the rates for each 
set of form factor, become higher when the central val- 
ues are compared: the difference provides us with hints 
on the level of improvement in the determination of the 
form factors needed at present. Such an improvement 
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FIG. 1: Differential branching fraction dBR(B s — » r]^£ + £ )/dq 2 {£ = e, /x) as a function of the momentum transfer q 2 (upper 
plots), and branching fraction BR(B a — > r)^£ + l~) versus the compatification scale (lower plots) obtained using set A and 
B of form factors. In each upper plot the lowest (red) region is obtained in SM (1/R — > oo), the intermediate (blue) region for 
the compactification scale 1 /R = 500 GeV and the uppermost (green) region for 1/R = 200 GeV. 



also concernes other determinations of the form factors 



As a further remark, we notice that the KLOE fit of the 



r] — rf mixing eq.(3.6) corresponds to using an effective 
mixing angle ip e ^ — (44.5 ± 1.5)° in the expressions of 
the B s — > rf form factors, and to a 10% decrease of the 
corresponding decay rates. 



ft is worth comparing the predictions ( 4.4 )-( 4.5 1 with 
the results quoted by the Heavy Flavour Averaging 
Group (HFAG) for the analogous B decay modes [36] : 
BR(B° -> K°e+£-) = (3.2tg;l) x 1CT 7 and BR(B+ -> 
K+l+l-) = (4.9 ± 0.5) x 10~ 7 for £=e,fj,. On the other 
hand, in the approach based on the QF mixing scheme, 
a relation connecting the decay width of B s — > r]£ + £~ to 
that of B s — > K°£ + £~ induced by the transition b — > d: 



1 



BR(B S -> K°£+£-) _ 
BR(B S -> r)t+l-) ~ (sin 



V td 



(1+A 



5(7(3)' 



-A') , 
(4.8) 

involves an 5/7(3) correction term Agjj^ (presumably 
small) and a term A' coming from the additional con- 
tribution proportional to V u bV* d , the correspondent of 
which has been neglected in the effective Hamiltonian 
(2.1 ) in case of b — > s modes. Such a kind of relations is 



of great interest, both from the theoretical and the ex- 
perimental viewpoint, and deserves a dedicated analysis. 

In the ACD model the Wilson coefficients C^* , Cg 
and Cio depend on the compactification scale l/R, 
hence decay rates and distributions vary with this pa- 
rameter, depends on 1/R through two functions 
D'(xt, jj) and E'(xt, jj) representing the contribution of 
7-penguins and chromomagnetic penguins, respectively. 
The contribution of Z°-penguins enters in the functions 
Y(xt, jj), Z(xt, ^) and X(xt, ^): Y and Z determine the 



coefficient Cg, while Cio depends only on Y. As for the 
function X, it enters in the ACD expression of the coeffi- 
cient Cl relevant for the modes with neutrinos discussed 
in the next Section. 



These functions obey the general representation (2.51, 
and can be found, e.g., in [51 [3] ■ For values of the com- 
pactification scale 1 /R of a few hundreds of GeV the coef- 
ficient Cy is suppressed in the ACD model with respect 
to the SM value, while Cio is enhanced and Cg turns out 
to be almost unaffected. For example, at = 300 GeV 



one has: 



c eff,ACD 



0.82 



C eff,SM 

1.16. As a consequence, for the d ecay into rj and mg = 



(jACD 



1.01, 



r<ACD 

r>SM 



we find for the function p(q ) in (4.2) at q 



1/R 



300 GeV- P( ^ aa 



|l/fl=300 GeV 



P{lliax)\sM 



1.23 - 1.42, de- 



pending on the set of form factors. 

The dependence of the rates on 1/R is obscured by 
the uncertainty of the hadronic form factors, as shown in 
Fig. [T] at values of this scale excluded by the analyses 
of the electroweak constraints: 1/R > 250 GeV. Sen- 
sitivity of the decay rate to 1/R could be achieved by 
a reduction in the hadronic uncertainty of about a fac- 
tor of two: in this case, the measurement of the rates 
of B s — > rj^'£ + £~ would allow to exclude the 1/R re- 
gion below 350 — 400 GeV, keeping the same uncertain- 
ties quoted for the other input parameters. Concern- 
ing the differential distributions ^{B s — ► rj£ + £^) and 

jj^i{B s —* rj'£ + £~) also depicted in Fig 1 the enhance- 
ment at small values of dilepton invariant mass obtained 
for decreasing compactification scales 1/R could be ob- 
served after the same improvement of the accuracy of the 
form factors. 
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A possibility to reduce the effects of the hadronic un- 
certainties consists in comparing the decay distributions 
of B s — > rj£ + £~ and D s — > i]£ + V( (£ = e, /i) at large val- 
ues of momentum transfer q 2 ~ i.e. close to the 
end-point. In this range the ratio: 



R(v) 



dT(B s 
dT(D s 



r)£+£-)/dy 
> r]£+vt)/dy 



(4.9) 



obtained evaluating the distributions at the same value 
as: 



can be written, up to corrections Of— M, 



R(v) 



v tb v t * 



V r , 



2 2 

a 



8tt 2 \M d 



\C 



10 



Co 



2(m b + m s ) e ffF T 



(y) 



(4.10) 



differential decay rate — — reads: 
dx 



dT(B s -> rfJ'vv) 
dx 

_ o \°l\ 2 \Fi s ^ ,l( {g 2 )\ 2 / 3 . 2 2 2 : . . 
" d 487r 3 M Bs V A (9 ' A7 b.' A V)J ' (5 ' lj 

with the invariant mass of the neutrino pair expressed 
in terms of x: q 2 = M% (2aj-l)+M 2 w . C L is the Wilson 



coefficient in (2.3 1 and the factor of 3 corresponds to the 



sum over the three neutrino flavours. 

The differential branching fraction is plotted in Fig. [2] 
in the case of 77 and rj , respectively, for the sets A and 
B of form factors. In the Standard Model one predicts: 



BR(B S -> i-jvv) = 



(0.95 ± 0.2) x 1CT 6 set A 
(2.2 ±0.7) x 10~ 6 setB 
(2.9 ± 1.5) x 10~ 6 setC 



(5.2) 



thus providing an access to the Wilson coefficients CV^ 

and C10, since the ratio of T Bs ^y j is practically inde- 

pendent of y and of the set of form factors. To further 
reduce the impact of the hadronic corrections, the double 
ratio could be considered: 



R(v) 



dT(B s -» r,£+£-)/dy 



dTiD^^+^/dy 
which is given by 



dT(B s -> K£+v)/dy 



dT(D s 



K£+v t )ldy 
(4.11) 



R(y) 



v tb v t : 



Vr... 



V, 



cd 



\c 



in 



Co 



2 2 



8tt 2 M 



M 



B.. 



2(m b + m s ) e ff F T S ^ V (y) 
M Bs +M v Cj F^ n (y) 



(4.12) 



up to corrections 0(m s ( =- )) E3- However, the 

\ Tfl c iii'b J 

double ratio involves many CKM suppressed transitions, 
and its measurement is challenging. 



BR{B S 



(0.9 ±0.2) x 10~ 6 
i] vv) = { (1.9 ±0.5) x 10~ 6 
(2.4 ± 1.3) x 10~ 6 



set A 
setB 
setC 



(5.3) 



These predictions can be compared to the present up- 
per bounds for the corresponding B — > Kvv modes, 
quoted by HFAG: BR(B° -► K°vv) < 160 x 10~ 6 and 
BR(B+ -> K+vv) < 14 x 10" 6 [36 , which are compati- 
ble with the SM expectations [3"B"] . 

In the ACD model, the coefficient is slightly enhanced: 

qACD 

for example J; SM ~ 1.10 at l/R = 300 GeV. How- 

Cl 

ever, the l/R dependence of the branching fractions is 
obscured by the hadronic uncertainty, as shown in Fig. 
[2] The distribution in missing energy shows an enhance- 
ment at small values of x which becomes sizeable at the 
lowest values of the l/R range, l/R - 250 GeV. The 
presence of two neutrinos in the final state makes the 
measurements of these modes a challenge that can be 
faced at an e + e~ superB factory operating at the T(5S) 
peak. 



VI. THE bs UT TRIANGLE IN ACD 



B s 



r\yv AND B s — > r/vv 



The exclusive decays with two neutrinos in the final 
state, which from the theoretical side are among the 
cleanest FCNC processes, require only one hadronic form 
factor. The observables are the decay distributions and 
the total decay widths. For the former, it is convenient 
E 

to use the variable x — — , with E miss the energy of 

Mb s 

the neutrino pair in the B s rest frame (missing energy), 



with 0.5 < x < 1 



In terms of the variable x the 



To complete our discussion about B s physics in the 
ACD model, we turn to the bs unitarity triangle (UT) 
and, in particular, to the weak phase j3 s defined as: 

V ts V T l~ 

Being ACD a Minimal Flavour 



Ps = Arg 



V CS V 



cb 



Violation model, the CKM matrix has the same struc- 
ture as in the SM: it is constrained to be unitary and 
is described in terms of four parameters, one of which is 
a complex phase: p, rj, A and A in the Wolfenstein pa- 
rameterization 39J . At 0(A 3 ) only V u b and Vtd have a 
complex phase, while at 0(A 4 ) also V ts is complex and 
its argument identifies (3 S j4"0] . 
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200 400 600 800 1000 200 400 600 800 1000 200 400 600 800 1000 200 400 600 800 1000 
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FIG. 2: Differential decay rates dBR(B a -> ^ r, w)/d X ELS £L function of x = E miaa /MB S (upper plots), and branching fractions 
BR(B a — » rf'^vv) versus J? -1 (lower plots) obtained using set A and B (right) of form factors. In the upper plots the lowest 
region is obtained in SM; the intermediate (blue) one for 1/R = 500 GeV, the uppermost (green) one for 1/7? = 200 GeV. 



As it is well known, the various unitarity conditions 
are represented as triangles in the complex plane. Even 
though the CKM matrix has the same structure in ACD 
and SM, and, in particular, its phase is in both cases the 
only source of CP violation, it may happen that the vari- 
ous CKM elements differ in the two models and therefore 
the UT triangles do not coincide. Deviations could be 
expected for elements extracted from loop-induced pro- 
cesses, where the tower of KK modes may play a role; for 
elements obtained from tree level decays no difference is 
expected, as for \V US \, \V u b\ and \ V c b\. For the db triangle 
defined by the relation: V ud V* b + V cd V* b + V td V t * b = 0, 

only the side Rt — J^ d ^* fc l depends on R^ 1 through 

\v cd v cb \^ 

the dependence of \Vt d \ coming from the KK contribu- 
tion to the — mixing [8 . Using such a result, 

together with the measured value of | " b J , the variation 

I Kb I 



of p — \1 — 4j- j p, r) = \1 — ^-J r], and of 7 (the phase 

of V* b ) has been obtained [H] , finding that the db triangle 
in the ACD model could be slightly different than in the 
SM. 

The bs triangle stems from the unitarity condition in- 
volving CKM elements relevant for B s decays, and is de- 
fined by the relation: 

V us V* b + V cs V; b + V ts V t * b = . (6.1) 

Analogously to Vt d , V* s inherits the R" 1 dependence from 
B° s - B° s mixing 0: AM S cx S(x t , l/R)\V ts \ 2 , where 
the proportionality factor does not depend on R^ 1 and 
is common to ACD and SM. The function S(xt , 1 / R) 
takes into account the SM contribution, So(x t ), and the 
contribution of the KK modes: S(x t ,l/R) = So(x t ) + 
J2n S n (x t , 1/R), with: 



c ( \ gt(5-10Jf + zf) irMwR f x t (Zx t -I) (1 - 3x t + lx\ - x\) 

\ Sn{xu Xn) = 8{ Xt - if + ~2~ 1 " %A - 1) 3 ( ' , ] ( ' 1 ] 

' + *" ' ~ } - J(R, 3/2) + -L \ dy [x t + 3(1 + x t )y]y/y ctaohfrMw Jtyy)} (6.2) 

4Xt Jo > 



2x t (x t - iy 



r 



and As a consequence one finds: 



-i v v t: ' 

J(Rv) = l d y yictanh( 7 rM w i?^)- a ;, 1+ "ctanh( 7 rm t i?^] i . F . gr 3j we ^ ^ dependence rf % 



ts\ACD on R 

(6.3) using the quoted value of |V t;s |sM- Since at 0(A 4 ) 
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FIG. 3: CKM element V ts (left) and phase f3 s (rad) (right) versus RT 1 . 



-AX 2 



Ps — arctan 



TCI' 
2A 2 t? 



2p 



2«ry], it turns out that 
Using the dependence of 



2- A 2 (l -2p) 
p and 77 (or p and 77), also the /3 S dependence on Br 1 
can be obtained; it is depicted in Fig. (3) The SM result 
for f3 s is small: (3 S ~ 0.019 rad, and the dependence on 
the compactification scale further reduces this value. 
As mentioned in the Introduction, preliminary Tevatron 
results point towards larger values of (3 S , so that our anal- 
ysis supports the conclusion that a confirmation of the 
measurement of a large phase in the B s mixing would 
point towards new Physics models different from MFV 
scenarios. 

In Fig. [3]the rescaled bs UT triangle in the ACD model 
is displayed at the representative compactification scale 
R- 1 = 300 GeV. 




0.2 0.4 0.S 0.8 1 1.2 1.4 



FIG. 4: Rescaled 6s UT triangle in the SM (continuous line) 
and the ACD model at R' 1 — 300 GeV (dashed line). 



The results for the branching fractions of modes with two 
charged leptons are O(10~ 7 ) for t = e, fj, and 0(1O -8 ) for 
£ = t, suggesting that they are within the reach of facil- 
ities like a SuperB factory. In the case of the final states 
with two neutrinos, the branching ratios are of O(10 -6 ). 
The uncertainty in the numerical results is dominated by 
the error of the hadronic form factors: in order to be sen- 
sitive to compactification scales 1 /R ~ 350 — 400 GeV in 
the ACD model the error in the form factors should be 
reduced by about a factor of two. 

A remark concerns the issue of the 77 — rf mixing. In 
addition to modes with light mesons, processes involv- 
ing heavy mesons provide information on this mixing. 
Examples are the J/ip radiative decays: J/ip — > 7/7 
and J/ip — » 777, which are sensitive to the glue content 
of 77 and 77' [41] . Although experimentally challenging 
and considered mainly for different purposes, also the ra- 

. BR(B S -> rt't+l-) 

tio „ -7 „ „ , „ , represents a possibility to dctcr- 

BR{B S -> ril+e-) F 

mine ip and to look at the gluonic content of 77'. A dis- 
crepancy in the comparison of the determination of the 
same quantity from other channels, e.g. from the ratio 
BR(4> — > 77 7) nj „ a -„„ i „ D , „(/) 



BR(4> 
modes. 



777) 



would signal peculiar effects in B s 



77V 
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VII. CONCLUSIONS 



Together with other rare B s decays, the loop-induced 
transitions B s — > rj^£ + £~ and B s — > rf'^vv must be in- 
cluded in the Physics programmes of experiments aimed 
at performing tests of the Standard Model and at search- 
ing signals of new Physics like those related to extra di- 
mensions. Within SM, using the flavour scheme for the 
77 — 77' mixing and the angle determined by the KLOE 
Collaboration from radiative <p decays, we have predicted 
decay widths and distributions of these rare B s channels. 
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